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Current Results for the Proposed Research 197^"75 

The investigation of the shock structure in an ideal dissociating 
gas based on the formulation of continuum gas dynamics has been carried 
out. Please see details in the attached paper, "Shock-Wave Structure 
in an Ideal Dissociating Gas." 


Proposed Further Investigations 197 5 “T ^ 


A, Investigation of "model dependency" on shock profile. 

For simplicity, Ligh thill's "ideal dissociating gas" model is used 
for the current investigation. It will be interesting to test the 
dependency of this model for the obtained results. 

B . Kinetic Theory Approach 

The current analysis of the shock-wave structure in an ideal dissociating 
gas is based on the continuum gas dynamics. It will be interesting to 
use the kinetic theory approach to investigate the same problem and then 
compare their results . 



Shock-Wave Structure in an Ideal Dissociating Gas* 

K. H. Liu 

Department of Physics, Southern University^! Baton Rouge, Louisiana 7^813 


The shock-structure problem utilizing Ligh thill's "ideal 
dissociating gas" model is investigated. The analysis is based 
on the usual formulation of continuum gas dynamics. The two 
extreme cases of chemical equilibrium and chemically frozen flow 
are investigated. Profiles for species concentration, temperature, 
and mean velocity are obtained for both limiting cases at various 
upstream conditions. The problem is also solved for arbitrary 
chemical reaction rates. Mathematical nature of the upstream 
and downstream singular points has been investigated and shock 
profiles for two typical chemical reaction rates are presented. 

As far as shock thickness is concerned, the solution for frozen 
flow shows good agreement with results obtained for a pure gas. 
Also, shock thickness based on the maximum slopes of temperature, 
species concentration and mean velocity, respectively, are 
practically the same; this is in accord with the conclusion of 
Talbot and Scala. 


*Uork performed in part, under the auspices of National Aeronautics and Space 
Administration under Contract No. NGR 19-005-003. 


Nomenclature 


velocity of sound 

specific heat at constant pressure 
specific heat at constant volume 
binary diffusion coefficient 
dissociation energy 
thermal diffusion coefficient 

'specific internal energy 
specific enthalpy 
forward reaction rate constant 

reverse reaction rate constant 

—LeL. , see equation ClO) 

RnT 

mean free path 
characteristic length 
Mach number 
pressure 

Universal gas constant 
temperature 

diffusion velocity of i^ species 
mean flow velocity 
molecular weight 
space coordinate 

"fctl 

mole fraction of i species 


Greek Symbols 


c<i = mass fraction of i-th. species ( 

Y = adiabatic index 

T - characteristic reaction time 

A = shock thickness 

k = compression viscosity 

/* = shear viscosity 

y ~ frequency of the vibration mode 

f 5 = density 

f 3 = characteristic dissociation density 

cb t . = ^net rate of production of i-th species 


Subscripts 


( )q conditions at either x = ±co 

( )i conditions at x = (upstream) 

( )2 conditions at x = +co (downstream) 

( ) e Equilibrium Flow 

( )f Frozen Flow 

( )^ i-th. Species 

i=a atomic species 

i=m molecular species 


over a large temperature range (lcA^xlO^K) , Lighthill asserts 
that is practically a constant. Such a gas is termed as an 
"ideal dissociating diatomic gas," its rate of production for the 
atomic species, to , is found to be 



with 

r = / [ 4k r p 2 0 

V can be interpreted as a characteristic reaction time. If the 
chemical reaction occurs very rapidly, the value of oi will differ 
'only slightly- from the local equilibrium value. In the limit of 

infinitely fast reaction, H « 0, the dynamical change of the flow 

* 

field is very slow in comparison with the dissociation and recom- 
bination rates. Since the reaction rate is reversible, *the flow will 
be isentropic^ . At the other limit of large V (very slow reaction), 
the rate of production, co , is effectively zero. The gas may then 
be considered as chemically inert and is usually referred to as 
chemically frozen flow. Again, the irreversibility associated with 
the reaction rates disappears and the flow, as in chemical equili- 

fn } 

brium is, in fact, isentropic w ' . 

In any analysis of chemical non-equilibrium flow, the main 
difficulties arise from the coupling between the gas dynamic equa- 
tions and the chemical relaxation equation. However, in the two 
extreme cases of frozen and equilibrium flow, the chemical relaxa- 
tion equation is reduced to a very simple form, thus greatly sim- 
plifying the analysis. For the non-equalibrium case, the assumption 


8 . 


of F being a constant will be used throughout the present work. 

(10) 

This assumption is consistent with the one made by Freeman 


Transport Properties 


The solution of the governing equations requires a knowledge 
of the variation of the transport coefficients 5 shear viscosity ^ , 
thermal conductivity X , binary .diffusion coefficient D am , and co- 
efficient of thermal diffusion D^., with temperature and concentra- 
tion;. For an actual engineering problem, experimental data may be 
used; but, for the case of an idea dissociating gas mixture of dia- 
tomic molecules, a simplified kinetic theory approximation can be 
made. 

Following Dorrance^^, we approximate the viscosity and 
thermal conductivity by 


^ = £ a ys . = Z KYi 

i <• 

with yj_, the mole fraction of each species. For a binary mixture 
of atoms and molecules, yu and X are 


t ~/ i * 

A 


( ~ oL 
i ■+ oL 

I - oL 

i + 


+ /«a 


z* 


i 4- ct 

ZcL 
| *4- eC 


respectively. Accurate to the first order approximation in solving 
the Boltzmann equation through the Chapman-Enskog procedure, we may 
further prescribe that 

. /*m -/‘a, 


JX -/*■ 


thus , 
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A first order approximation to the thermal conductivity 

( ip} 

for a monotomic gas is given by Penner' * 

^ &/«■ 

For a pure polytropic gas, the thermal conductivity may he appioxi- 


mated 


( 12 ) 


by 


v - u 4 ^ 

Aw “ A K »/ t » v. 3 15 R n y 


The factor in parenthosis is known as the Eucken approximation which 
accounts for the transfer of energy between the translational and 
rotational degrees of freedom when polytropic molecules collide. 
Approximating the specific heats at constant pressure, namely 
C^ m and C^ a , by their values for chemically frozen flow, one finds 

from the caloric equation of state that 

♦ * 

= otC^ + otnC^ - [ K„ + (5~ K»)*] Rn ' (16) 

hence 

X = -4^- ■ [ (5 + 4 K„) + C55- 4 K")" 4 -] 

4O +c 0 v 

x C5+4KO + (55- AK n y 
= C17) 

Eq.(l7) implies that, although the mixture viscosity fi may be 
assumed the same as that for the atomic species, the thermal con- 
ductivity X will be less than X«. by a factor of about 2/3, since 

* l ■' 

Kj^ varies between 7/2 and 9/2. ’’ - 

The Prandtl number Pr , through combining Eq,(l6) and (17) , 


becomes 


10 . 



4(f + «0 f ^n + C5- 

£ 5” + 4- K w ^)t 4 5?" 4 14/)^ 


( 18 ) 


Taking the average value of K m , Km— 4, one finds that for oL vary- 
ing from zero to one, Pj* from Hq. ( 1 8 ) varies from .16/21 to 40/60. 

Using the thermal equation of state to express density in 
terms of other thermodynamic variables, we have the Schmidt number 
Sc . 


Sc 


p r^<3w 


f 1 Oam 


( 1 4 °Q Rn) T 


Assuming that fDo.-n is independent of composition, it may be eval- 
uated at cL-0 , Hence, we obtain 


Sl = (-pk~) ' 0 + 0 “ O«0 ' <19) 

V Pm iV ' 


The thermal diffusion coefficient D-j. may be expressed in 
terms of the binary diffusion coefficient V m . 




(20) 


Theoretical results for kt based on kinetic theory are very complex 
but to an acceptable approximation 


k+ - 0. 3 5 


Wm - V/g 

VU + YL 


)4 


In the case of a pure diatomic gas 


thus from Eg. (20), 


kn 


0.234 


c> - *y 


Dt 


0.23+ 


ot ( 1 - «Q 
(t-f erf.)* 


I^om 
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The coefficient D^/pD 0fn appearing in the diffusion velocity [see 

Hq. (24)] may be written as - 
* 

D t = ■ W a W„ , 

r * 


Since 


therefore 


WgWm _ Q + oQ ; 

iT/ 4 


~^r- - 0.117^0-^). 

Governino Equation in Dimensionless Form 


Governing equations are normalized by introducing the fol- 
lowing set of dimensionless quantities: 


V = wu/p 

(velocity) 

T = TO 1 R„T / P* 

(temperature) 

A = wQ/P z 

(shock strength) 

W*D,n/P l 

(dissociation energy) 

f - m r/L W w 

(chemical reaction time) 

X = A + -£- 
.3 

(viscosity) 


and the characteristic length, L s JJ. / 'rt\\ 

Thus, the continuity equation for the atomic species, the momentum 
equation, and the energy equation become ^ 

+ [O-O-jfe"'*- * X J (an 




(21) 

( 22 ) 


[ (K n - -!• " ^ j » v V J> 


* • > t 

tj — t + w 


1 f + (5" Jt 

„ . t— W ** P .. . . . 


" d* 


+ [C$-Km)^ + Tj>]Jk^ = 


12 . 


respectively. In Eq.(23) the term due to the diffusion velocity 
(a^/m is given by 


M*- ( 1 ' \ f r | , ( 1-^) 1 , r_Pt *0-0 ] La-g 

“ \X/p OR™ x. > 


X dX 


4* 


- ccQ-^'L d_v 

Z v dx. 

This set of equations for v, % and oi ] the boundary conditions pre- 


(24) 


scribed at upstream and downstream of the shock v/ave 

dV _ dT dri 
dX 


-0 


at X «± 


03 


dX dX - dX 

and the transport properties of the gas: \ f ft , II., Dam, X completely 
define the shock structure problem. 

The set of governing equations for the three cases that we 
shall investigate can now be listed as follows: 

(A) Equilibrium flow (F = 0) 

Species Continuity Equation : 

Since F =0 in Eq.(21) the term in brackets must also be 


zero, thus .the equilibrium mass fraction is given by 

.2. p _ IS 

* „ Vp - -v 

I - ot " P 

which is merely the repetition of Eq.(l4). 


( 2 ?) 


Momentum Equation: 


dv 


XL ^T " v + 0 + V 


1 


(26) 


Energy Equation : 

Substitution of Eqs.(24)(25) into (23) gives 
{ £ t K„ + (5-K„y) + 


V* 


~ C Kri~ 0 t ( 4 - \ -!■ * ' 


£1 
/ 1 




( 27 ) 
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(B) Frozen Flow ( F -»-«0 

Species Continuity Equation: 

L dv 


lf — — m - t C «-«o 

_ r P t , *'o-O i _l 

L P D ^ J V 


‘V"2T + 


f Am 

Momentum Equation: 
dv 


2 

ri fl-sp l _L dv_ 
dx 


XLg = v+ 


(28) 


( 29 ) 


Energy Equation: 

J r + JL^| = [ (f-Km)*!- ^Jt + V + 


v z . ^ T 

“ y t T J) 2 


( 30 ) 


(C) Non-Equilibrium Flow 

The terms related to the diffusion velocity in the 
species continuity equation and the energy equation will be neglected 
in the case of finite chemical reaction time. The simplification due 
to this assumption is obvious, since the species continuity equation 
is hence reduced to a first order differential equation. The assump- 
tion is justified from the results for the frozen flow case where the 
chsnge in concentration, caused only by diffusion, is found to be 
very small. In addition, the chemical reaction time T is taken to 
be constant throughout the shock wave. Then, the set of governing 
equations becomes 


l = -L [ e 

u dx r 1 v 


f* 


( 31 ) 


XL jj~ = v + (1**0 “ 1 


( 32 ) 


lb. 


-£* [ Km + C 5*- Km> } t- 55 f( K ^‘ 0 4 'P+eiT,-^.+\/-A (33 ) 

Mach Number 

The basic parameter in the shock Structure Problem is al- 
ways the upstream Mach number, M = u/a, in which the speed of sound a 

is ^rJu^OKJ 

Since the adiabatic index, y , is not defined in a simple way if con- 
centration changes are allowed. 

For the case of frozen flow 

'■ ' Km + (5- ftQot 

^ (Km- t) + (4-Km)U 

In the. case of equilibrium flow, Ciarke^^) gives the ratio Y<®/fr for 

d(f-U) ( ('Vl*00+«0“ (4 +*Q] Z ) " { 

3 (Tp+ i/d (\-d x )+2(4*°0 J 

One can see then 

r e ^ r f 

Since the frozen and equilibrium situations represent the two extremes 
of chemical behavior, the adiabatic index for the non-equilibrium flow 
will be bounded between and Y e * We have, therefore, three different 
speeds of sound making the standard definition of the Mach number some- 
what ambiguous. To circumvent the complications that could arise from 

x 

this ambiguity, we adopt the flo w parameter YM which is applicable 
to all three flow situations. 

z v ~ v 

y — y 

~~ (1+ et )R„T ~ 


to be 



At upstream as 


X — > - CO, 


(VM l ), = 


'I 


1 0+°<i)Ti 

Substitution of the momentum equation gives 




I- v, 


or 


v l = 


(YM 1 -), 

1 + (y 


Therefore, iii the present analysis , or V, will be used as the 

basic parameter. Herein, this parameter is assumed to lie in a 
range where the existence of a compressive shock wave is assured. In 
other words, the strength of the shock wave is large enough to give a 
continuous shock profile for the all values of P . 


Rankine-Huaoniot Relations 

Since the- gradients of all physical quantities vanish at a 
large distance upstream and downstream of the shock wave, conditions 
at end points for equilibrium and non-equilibrium flows are the same. 
The end consitions for frozen flow are different due to the fact that 
the concentration downstream of the shock wave remains the same as 
the upstream value. 

(A) Equilibrium Flow 

Setting all gradients zero in Eqs. (26,27), we find 
the norma] ized velocity at end points to be 



I + pc 



(34) 


O — 


c e 


K„ f (5- 


CA-I'A (a 


— } K. 


(■ 


> 4 .. 

“ I % *1^ -V 
1 ^ 


__ 2 (. V& ~ ' J f _ 


A.* 


^ z 


i 


where 
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The plus sign in Eq. (34) is associated v/ith the upstream point, 
while the minus sign should be taken for the downstream point, x«-po. 
Through the momentum equation, Eq.(26), we find the normalized temper- 
ature at end points v to be 

_ _ ^ C 1 ~ v jiil 

^ _ — ; “ 

I + (35). 

The concentration oi^ is then found by knowing \ and %, through 

the species continuity equation, Eq . (25) • Eqs. (25434,35) are solved 

simultaneously to determine v, X , and ci a t the two state points. For 

illustration. Figs. 1,2,3 present the results for shock waves of 

various strength characterized by (VM*),. It is interesting to note 

that both the velocity ratio and the temperature ratio for a given f 

upstream Mach number do not vary monotonically with respect to the 

* 

change of upstream atomic concentration. 

' (B) Frozen Flow 

; The conditions for the frozen flow case differ in 

I 

j that the mass fraction o4 is the same at x =i«> . Since the upstream 

■ conditions represent an equilibrium state, the downstream state 

I 

j conditions ( do not in general represent an equilibrium 

. state. Thus, the flow after passing through the shock must make a 

' transition to an equilibrium state behind the shock. This relaxa- 

! tion zone, analogous to shock transition in a. two-phase (solid-gas) 

i fluid, will not be discussed in the present analysis. Readers are j 

(tuv 

referred to^text book on physical gasdynamics for further discussions. > 

* j.< 

lj 

Hence, instead of having Eq. (2?) for the concentration, we now have i 

the condition 

’ O^j — Ci i. r — C*.q I 



17 . 


together with the momentum equation and the energy equation. Setting 
all gradients to zero, we obtain the normalized velocity at both end 
points 



(36) 


where 


P f = 


Km » (£_- KiQ^o 




Again, the plus sign should be adopted for x = -«> ; the minus sign 
for x = +<* . The normalized temperature at two end points is 

T = Cf" Vi,z ) 


I A 


I + oL t 


A comparison of Hq. (3*0 and Hq. ( 36 ) shows the only difference in V itZ 

between equilibrium and frozen flow is the replacement of in the 

former case by in the latter case. Since the concentration oL 

is unchanged for frozen flow, is thus constant, the values of all 

physical variables at end points for frozen flow can therefore be 

obtained without numerical iteration. The parameter, E* , which is 

also constant may be written as 1 

^ [ M fi - 


1 


f ~ . 1 i (4- rMf , J 

in which M* ( is the upstream frozen flow Mach number. 


Here 

and 


= 0 for M’j. j = 1, 


c-~ — — • 
' j, _ 


1 or M 


i 1 " 


°' <E ^7TT 


thus 
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Phase-Space Analysis . 


It is a characteristic feature of the shock structure 
problem that state points at upstream and downstream are unavoid- 
ably singular points with regard to the set of governing equations. 
Phase-space analysis is the usual technique applied to assure the 
uniqueness of the solution mathematically, and in practice to ob- 
tain an initial slope to start the numerical integration. Analysis 
in the two dimensional phase plane has been a standard procedure, 
while the investigation in three or more dimensional phase space j.s 
by no means a simple routine. The classical treatise of Poincare^^ 
defined the nature of singularities in three dimensions, but did not 
illustrate the procedure with practical problems. The present 

analysis follows the path ’set forth by Talbot and Sherman^ in 

* 

their investigation of the shock structure problem based on continuum 
gasdynamics, taking into account the higher-order Burnett stresses. 

For an autonomous system of three non-linear ordinary dif- 
ferential equation governing v, X , and , we have 

A dT = 

B n = s c v, -o) 

c ai = H (S*'. 



where A, B, C are simply constants, which, if one prefers, may be 
absolved in F, G, and H, respectively. The boundary conditions at 
end points are such that 

F V S/ T P ) » G % *&») - H (**, V*, ^0 = 

thus making the slopes , 4 / 4 ■/ indeterminate at end points. 
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Introduce 


£* ss oi — C < 0 ' 


Vr 


v * V - Va *» v 0 <2. 


her 


^> 3 - 


J 


(38) 


X - -X - % ~ & 

fs/ ^ 

with «£>, Vo,^ being constants, (T a parameter along the integral 
curve, b a characteristic root. In the neighborhood of the singular 
points, Eqs.(37) and (38) give us a system of algebraic equations 

A/ ^ ^ -i 

governing ot 0; . Vi T> . 

- Ab -J* Fa 

Bb+ Gt/ 




-Cb+Hi 


* 



/N/ 

Vp 


A/ 


T *J 


- w 


(39) 


He* H v 

where = (|§) o f FV = (pv") » etc * T° have a °f non-trivial c? 0 ,% 

* 

and T, , the determinant of the square matrix in Eq. ( 39 ) must 
vanish, thus resulting in a cubic equation for b. 

(ABC) b 3 - [ACGy + AB H v + BC Rj + 

+ f A (<5v H»r H*> B (F* - F p fU + C 05ftr b+ (W 

+ ( Fv(<5o(Hf “GtH^) + fu (G-t Hy-§vHt) + F-c - G 

In general, three different roots of b are obtained, among which one 
value of b should be selected; according to which a set of d 0) % X 
is computed. The selection of this particular b is illustrated in 
the following calculation. 

(A) Equilibrium flow 

Since in this limiting case the mass fraction may be written in terms 
of the velocity and temperature, the number of independent variables 
is reduced to two. Thus, in the set of equations [Eq.( 3 ?)j v/e have, 
equivalently, 


i 


F 0*,v,tO=o. 


Or <L.&S = H C‘tl) 

3 dv G 

One can also shew that the cubic equation, Eq.(4o), is reduced to a 
quadratic. An analysis of the two roots, bj, b 2 , of this quadratic 
at each of the state points reveals that the downstream point, , 

is invariably a saddle point; and the upstream singular point, x’-s° , 
is an unstable node. In other words, the ratio b i/bz is 

negative, thus signifying a saddle point; it is positive at , 

thus signifying a nodel point. 

Take a typical example where the upstream atomi* concentra- 
tion c<, is 0.183; we sketch the variation of b as a function of the 
normalized upstream, or downstream velocity in Fig. 4. There is a 
downstream velocity v^ associated with every given upstream velocity 
v^. One root remains positive for all values of v^ ana V 2 , the 
second root passes through zero at a certain value of v where v^=V 2 , 
or the upstream Mach number based on equilibrium sound speed is ap- 
proximately unity. The second root is obviously the characteristic 
root of interest, since there should be no shock wave of finite 
strength at equilibrium Mach number of unity. One also observes 
from Fig»4 that the value of b reaches a maximum and a minimum when 
Vj >0.9* This behavior need not concern us, as the corresponding 
Mach number for V, >0.9 is so high that the model of ideal dissociating 
gas breaks down long before the situation arises. 

Having determined the state conditions at end points and the 
initial slcne. Eq. (4l) car, be integrated numerical. ' in 2 routine 


manner. 
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(B) Frozen Flow 

Unlike the preceding example where the algebraic equation 
(Eq.(h-O)) for b is reduced to a c.uadr tic, the present case has a 
cubic equation governing b. The b versus V, variation is presented 
in Fig.^ for two different upstream conditions; =1 'and c^=0. 

Three real roots of b are obtained, among them one is chosen based 
on the same reasoning as given in the preceding section. The selec- 
ted value of b used in numerical integration cuts across the b=0 
line at two different places for =1 and c*,=0; because the defini- 
tion of upstream Mach number depends on the upstream atomic concen- 
tration. Once a proper value of b is selected, the integration pro- 
ceeds from the downstream point toward the upstream and point. 

(C) Non-Equilibrium Flow 

Eq, (40) is again solved for two examples where the reaction 
time f is J and 2 respectively. Results are presented in Fig. 6. 
Three real roots of b are obtained, the selected value of b for in- 
tegration is the one which crosses b=0 at Mf t ^ 1. The integration 
starts from downstream end point and proceeds toward the upstream 
end point. 


Numerical Results 

(A) Shock Profiles 

Shock profiles for an example of C YM^^IO and <*,=0.183 sre 
presented in Fig.7j in which a very pronounced overshoot in concen- 
tration ratio is recorded, while the temperature profile varies only 
monotonically . As the reaction time T increases from zero (equili- 
brium flow) to a value of 2 and then 5; we find in Fic.o and 9 a 
temperature overshoot which increases with f . In the limiting 
case of frozen flow, T— >co , the shock profile again records a 
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monotonic variation for the temperature, but a high peak for the 
profiles of the concentration ratio. In this case, the concentration 
ratio changes purely due to mass diffusion, not because of chemical 
reaction as in ‘the case of f =0. 

(B) Shock Wave Thickness 

The most commonly employed definition of the shock thickness 
is 


where |dxi max is the maximum slope of the velocity profile and ) 1 5 2 
is the total change in v. For comparison with previous investiga- 
tors, it is convenient to obtain shock thicknesses in terms of mean 

free path evaluated at the upstream state point. The mean free path 
• » • , 
for Maxwell molecules is 




6v 


( 42 ) 


l - li. A) 

1 5 P, ,/ ztt ft T 


Assuming the temperature dependence of viscosity is given by* 

(if ay 

the ratio L/l ( is determined to be 

l . mi ji/il 


i. 


16 


(™^i 


M Vi 


where the characteristic length, L, as defined previously is L= J 


h 




■A 


Though S should 
Ms .'.veil mo i ecu! 
to account for 


be unity in consistency with the 
es, it is not uncommon to assign 3 
~he departure of real gas from the 


assumption of 
seme other vs 
Maxwell model 


iue 
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The expression for the shock thickness *n units cf mean free path 
becomes 

A _ ; 

i, " I clv__ 

' ifBSjT (T/r,y | 4(^) 

in the present analysis the exponent S was token to have the value 

O.p. 

Figure 11 shows the variation of the ratio ^ /& with (7AI ), 
for both chemically frozen and equilibrium flow. It is apparent 
that the value of for the equilibrium case is greater than that 

for the frozen flow condition at a specific value of ( . In 

other words, for a given upstream condition the shock thickness 
based on maximum slope is greater for the frozen flow case than for 
the equilibrium situation. For the equilibrium case, the thickness, A 
for a given value of (VM 4 ) first decreases then increases as the 
atomic concentration is increased. A similar type of behavior was 
previously noted in the determination of the upstream and downstream 
conditions in equilibrium flow. This suggests that for a given 
value of (TM X ) ( the thickness, A , has a maximum as the atomic mass 
fraction is increased from zero to unity. Fig. 12 shows that this is 
indeed the case for flow in chemical equilibrium. 

Fig. 13 gives a comparison of the results from the present 
study with previous investigators. As expected, the solution for 
frozen flow shows good agreement with results obtained for a pure 



gas 
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Talbot and Sea 



define a shock thickness based on any flow 


property, f. 


A f " 


^viax- £ 


wm. 


_df_ 




Fig.l4 shows the variation with (YM;^ of the maximum slope thick- > 
ness based on v,*V, and oi for both frozen and equilibrium flow. 

It should be noted that the definition of A^ in Fig,l4 is slightly 
different from the preceding equation, as the factor f, tV3X — f M j„ t is 
replaced by j ’Ft — f x j ; yet the quantitative discrepancy in based 
on any physical quantity, even if it has an overshoot in the shock 
layer, is always negligible. It can be seen that different proper- 
ties predict shock thicknesses of the same order of magnitude; 
this is in agreement with the conclusions of Talbot and Scala. 
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FIGURE CAPTIONS 
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Fig 1. VELOCITY RATIO *yr FOR EQUILIBRIUM FLOW 

Fig 2. TEMPERATURE RATIO ^ FOR EQUILIBRIUM FLOW 

Fig 3. CONCENTRATION RATIO ^ FOR EQUILIBRIUM FLOW 

Fig 4. CHARACTERISTIC ROOTS (bi, bg) FOR EQUILIBRIUM FLOW 

Fig 5. CHARACTERISTIC ROOTS (bi, b 2 , b 3 ) FOR FROZEN FLOW 

Fig 6.’ CHARACTERISTIC ROOTS (bi, b 2 , b 3 ) FOR NON -EQUILIBRIUM FLOW 

Fig 7. SHOCK PROFILES FOR EQUILIBRIUM FLOW 

Fig 8. SHOCK WAVE PROFILES FOR NON-EQUILIBRIUM FLOW 

Fig 9. SHOCK WAVE PROFILES FOR NON-EQUILIBRIUM FLOW 

Fig 10. SHOCK WAVE PROFILES F^R FROZEN FLOW 

Fig 11. SHOCK THICKNESS FOR EQUILIBRIUM AND FROZEN FLOW 

Fig 12. SHOCK THICKNESS FOR EQUILIBRIUM FLOW 

Fig 13. COMPARISON OF SHOCK THICKNESS WITH PREVIOUS INVESTIGATORS 
Fig l4. MAXIMUM SLOPE THICKNESS BASED ON d , T AMD V. 
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